Spin and orbital Hall effects for diffracting optical beams in gradient-index media 
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We examine the evolution of paraxial beams carrying intrinsic spin and orbital angular momenta 
(AM) in gradient-index media. A parabolic-type equation is derived which describes the beam 
diffraction in curvilinear coordinates accompanying the central ray. The center of gravity of the 
beam experiences transverse AM-dependent deflections - the spin and orbital Hall effects. The 
spin Hall effect generates a transverse translation of the beam as a whole, in precise agreement 
with recent geometrical optics predictions. At the same time, the orbital Hall effect is significantly 
affected by the diffraction in the inhomogeneous medium and is accompanied by changes in the 
intrinsic orbital AM and deformations of the beam. 
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I. INTRODUCTION 

Spin-dependent transverse transport - the spin-Hall ef- 
fect (SHE) - of classical waves and quantum particles 
is currently attracting growing attention in condensed- 
matter [H, high-energy 3, and optical [1, d i, d, 0, Hi 
physics. This effect appears under bending of the wave 
trajectory in an external potential and is closely related 
to such fundamental phenomena as the Berry phase, con- 
servation of the total angular momentum (AM) of the 
wave, and spin-orbit interaction. 

The optical SHE deals with evolution of Gaussian- 
type wave beams bearing intrinsic spin AM. A general- 
ization of this effect for higher-order Laguerre-Gaussian- 
type beams with phase singularities (vortices) has been 
put forward recently 



Hi 



As such beams carry 
well-defined intrinsic orbital AM [f^, an orbital-Hall ef- 
fect (OHE) appears under the bending of their trajecto- 
ries. 

Extensive studies over the past several years mostly 
considered the semiclassical trajectory equations tracing 
evolution of the center of gravity of a wave beam rather 
than the propagation of real extended beams. At the 
same time, the typical transverse shift of the beam's cen- 
ter of gravity due to the SHE or OHE is proportional 
to the wavelength and is rather small as compared to 
the characteristic scale of the beam deformations in an 
inhomogeneous medium. Therefore, it is important to 
give a picture of the evolution of realistic beams in a 
gradient-index medium, which includes SHE, OHE, and 
the diffraction processes. 

Below we provide such a description of paraxial opti- 
cal beams carrying spin and orbital AM and evolving in 
a smooth gradient-index medium. Although we consider 
Maxwell equations, our analysis can readily be extended 
to quantum wave equations describing evolution of quan- 
tum particles in external potentials. In particular, the 
OHE arises from the Laplace operator in curvilinear co- 
ordinates and is universal for any beams with vortices. 



II. MAXWELL EQUATIONS IN THE 
RAY-ACCOMPANYING COORDINATE FRAME 

Maxwell equations for the monochromatic electric field 
E in a gradient-index dielectric medium read 



E 



k^'^V (V • E) = 0, 



(1) 



where kf) ~ lj / c {lj is the wave frequency) and e = e (r) 
is the dielectric constant of the medium. Classical geo- 
metrical optics (GO) shows that in the short-wavelength 
limit the wave propagates as a classical point particle 
moving along the ray trajectory Yc = Vc (s) given by [isj 



t , t 
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(2) 



Here the overdot stands for the derivative with respect 
to the parameter s, which is the trajectory arc length, 
t is the unit vector tangent to the trajectory, V± = 
V — t (t • V) is the gradient in the plane orthogonal to 
the ray, and the subscript "c" means that the function 
is taken on the trajectory, i.e., Vj_£c = (^±£)lr=r ' ^t*^- 
Equations ^ define the central reference ray and realis- 
tic beams evolve in the vicinity of it. 

The wave-beam propagation is described using parax- 
ial approximation in the vicinity of the GO ray This 
implies the smallness of the two parameters: 



A 



Ail 



< 1 and ^2 



w 
L 



« 1, 



(3) 



where A is the wavelength, w is the characteristic beam 
width, and L ~ | Ve/ep"'^ is the characteristic scale of the 
medium inhomogeneity. We aim to describe the evolution 
of a paraxial beam keeping the terms up to the fi'^ order 
with respect to the combined parameter jj, — max (/ii , ^2) 
in Maxwell equations. Previousl y, t his problem has been 
solved with an accuracy of /ji^ [iJl, i-c, in the lowest- 
order approximation that describes diffraction but does 
not account for the Hall effects. 



2 




FIG. 1: (Color online.) Geometry of the wave propagation 
along a curved geometrical-optics ray. Depicted are: a ray 
accompanying coordinate frame attached to the unit vectors 
(ei,e2,t), direction of the inhomogeneity bending the ray, 
VxEc II t, and directions of the Hall effects orthogonal both 
to t and t (the large double arrow). 



Equation ([8]) can be projected onto the plane (Cij^) or- 
thogonal to the ray. In so doing, we notice that [16|| 
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■Eni, S||t = 2(Vi-Ej 



where we used Eqs. ([2]), ([6]) and ([7]). This yields the wave 
equation for the transverse electric field E^: 



(9) 
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h ds \ h ds 



l_d_ / 9E_L 



eEi 



X V± X E_L = 0. (10) 



In the vicinity of the GO ray, we introduce a ray- 
accompanying coordinate frame (^i, ^2, s) attached to the 
unit vectors (ei,e2,t). These vectors evolve along the 
ray: t = t (s), = (s), j = 1, 2, see Fig. [1] To make 
the (^15^2)5) frame orthogonal, (61,62) must obey the 
parallel transport equation along the ray [3, [l3| : 



(e. • t) t. 



(4) 



The Lame coefficients of the coordinates (f 1 ,^2,5) are 
equal to [13, Hlj 
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where ^ is the radius- vector in the (^1,1^2) plane. 

The electromagnetic wave is near-transverse in the ray 
coordinates: 



E — E 



S||t 



E,e, 



E, 



^^\E\ 



(6) 



From the equation V • (eE) — 0, stemming from Eq. ([T]), 
it follows that in the lowest-order approximation in /.i. 



V E 



Ei and £'11 ~ jfcr^ V± • Ej 



(7) 



where kc = k^^/el is the central wave number. 

Using Eq. jS]) and the first of Eqs. jT]), MaxweU Eq. ^ 
in the (^1,^2, s) coordinates takes the form 
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eE 
= 0. 



(8) 



As we will see, the last term in Eq. ([TU)) . which is of 
the order of /i'^, describes the SHE of light. This term 
originates from the combination of the polarization term 
V_L (V • E) and the Coriolis term 2ikcEui 



III. PARABOLIC-TYPE EQUATION 

Full wave Eq. (fTU)) can further be simplified to a 
parabolic- type equation via the WKB ansatz: 

E^ (s, = e (s) £-1/4 (s) W {s, C) e**(^), (11) 



where $(s) =^ kg J ec{s')ds' is the GO phase along the 


ray, 6 (s) ~ aei (s) -t- /?e2 (s) is the unit polarization vec- 
tor, 6* -e = 1, and W (s, ^) is the unknown slowly varying 
envelope of the wave. The components of the polarization 
vector is unchanged along the ray, a, /3 = const, which 
signifies the parallel-transport law for the wave electric 



field, related to the Berry phase |14l . Il7l |. Owing to 
Eq. (Ill), derivatives 6; || t do not contribute to the equa- 
tion Uni) for E_L. Note that only the last term in Eq. (fTU]) 
involves the wave polarization, and this term is diagonal- 
ized in the basis of circular polarizations: e = 61 + icre2, 
where ct = ±1 is the wave helicity due to the spin. Substi- 
tuting Eq. (fTTj) with a circular polarization into Eq. (fTO|) 
and retaining only terms up to the /i'^ order, we arrive at 
the parabolic- type equation for W: 



2ikc^ + AxW + kl 
ds 



W ^i(7 



-2ik^ i 



d Vx£. 

ds 2£c 



X ViM^ -t 



w-\k 



Vi_W 



[i-V^Ye, W. (12) 
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and in derivation of 



Eq. (fT2|) we used the Taylor expansion for e (r) 



Here Vx = d/d^, = | 



1.. 



V 



J Sc. 



The first two terms of Eq. (|T2|) represent usual 
parabolic equation, whereas the next terms in square 
brackets describe the influence of the medium inhomo- 
geneity on diffraction. All the terms in the left-hand 
side of Eq. ([12]) are of the order of they have been 
obtained in Ref. [li]. The right-hand side of Eq. ^ 
represent corrections of the /z^ order. 

As we will see, the first two terms in the right-hand side 
of Eq. (|12p . which are proportional to 'V±W, describe the 



SHE and OHE. (Note that the second term originates 
from the term in the Laplace operator and has a 

universal form for waves of any nature [isj.l These two 
terms can be eliminated by a simple transformation: 



(13) 



Substituting Eq. ^ into Eq. ([12]), we notice that ^ = 



9s 

— S ■'V \W and the above-mentioned terms are can- 
as 

celed when 



o = — I — X t 



i V±ec 



As a result, the parabolic equation acquires the form 



(14) 
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ds 2ec 
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(15) 



IV. SPIN AND ORBITAL HALL EFFECTS 

Equations p^ - (|15p are the central results of the pa- 
per. While Eqs. ([T^ and describe the transverse 
deformations of the beam due to the SHE and OHE, 
Eq. (11511 describe all other deformations caused by the 
diffraction in the medium. The cr-dependent term in 
Eq. (I14p is responsible for SHE; it takes the form as pre- 
dicted by the modified GO theories [1, i, H, @, [1] ■ Equa- 
tions lfT3|) and (|15p (which is independent of polarization) 
imply that the SHE produces a perfect translation of the 
whole beam in the transverse direction, without any other 
polarization-dependent distortions, see Fig. [21 

The OHE is more intricate and is described by the 
imaginary term in Eq. (|14p . To show this, let us first 
consider an unperturbed beam carrying a well-defined 
intrinsic orbital AM. Its field contains an optical vortex^ 
which is described by the structure 



Wi^ +isign(/)e2]''' =pl'lexp(i/(^) 



(16) 



where / = 0, ±1, ±2, ... is the vorticity, characterizing the 
value of intrinsic orbital AM per photon, whereas (p, ip) 
are the polar coordinates in the ^ plane. It is easy to 
see that small imaginary shift along, say, the axis: 
^1 ^ ^1 — ix, deforms the intensity distribution in the 
vortex along the orthogonal ^2 axis: 



oc 



e? + (6-sign(0x)^ 



P 



1 



2M^. 



shifted along ^2 on the distance sign I x, while the center 
of gravity of the vortex is shifted along the ^2 axis on 
distance — i.e., in the opposite direction, see Fig. [5] 
and cf. Ref. [lT|. Taking this into account, one can 
derive from Eq. (|14p the differential equation describing 
the shift of the center of gravity of the vortex ^B)) : 



^-'> = ^ (t X t) 



(17) 



From this equation follows that the nodal point = is 



where we substituted — t from Eq. ([2]). Equa- 

tion ([T7|) represents correction to the ray Eqs. de- 
termining motion of the center of gravity of a wave 
carrying well-defined spin and orbital AM. It is in 
agreement with the geometrical-optics predictions of the 
SHE and OHE |a], which are directly related 

to the conservation of the total AM in the problem 0, 0] . 

However, the /-dependent term in Eq. ()17|) is valid only 
until one can neglect perturbations in the beam shape 
caused by the diffraction in a gradient-index medium. In 
the lowest-order approximation, the diffraction-induced 
deformations are described by Eq. I|15p with the right- 
hand-side terms neglected [IJ]. Typically, the beam 
acquires elliptical deformations at distances comparable 
with the characteristic inhomogeneity scale [l^. These 
deformations do not affect the SHE, but they do af- 
fect the OHE, because elliptical deformations of a vor- 
tex beam dramaticall y c hange the intrinsic orbital AM 
carried by the beam [10| and the shift of the center of 
gravity of the vortex. In contrast to spin AM, the in- 
trinsic orbital AM of the beam is not conserved upon 
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FIG. 2: (Color online.) Schematic picture of the SHE and 
OHE transverse deformations of the beam with respect to the 
classical GO ray, Eq. ([2]). Shown are the cases of of circularly 
polarized Gaussian beams with cr = ±1 and I — (upper 
panels) and of vortex Laguerre-Gaussian beams with / = ±1 
and a = (lower panels). The intensity distributions are 
plotted in arbitrary units using transformation (|13p with some 
S ocaixt + ii, Eq. ([14)1 . 

the diffraction in a gradient-index medium. Indeed, the 
operator of the orbital AM, L oc —i^ x V^, does not 
commute with the operator in the square brackets in the 
left-hand side of Eq. 

Equation ([T5|) cannot be solved analytically in the 
generic case even with the neglected right-hand side. 
Therefore, it is impossible to determine analytically the 
OHE shift of a diffracting beam in a gradient-index 
medium. However, Eq. (jlSp can be integrated numer- 
ically in each particular problem. Then, according to 
Eq. (fT4|) . the OHE shift and deformation of the beam 
can be taken into account by introducing an imaginary 
shift ^ — > ^ — ik~^t ds at each step ds. It should be 
noticed that the diffraction of the vortex beam also de- 
pends on the absolute value of the vortex charge, \l\, but, 
in contrast to the OHE, Eq. (fT5|) is independent on the 
sign of the vortex, sign(/). 

To conclude, we have derived a parabolic-type equation 
which describes propagation and diffraction of paraxial 
beams in a gradient-index medium and accounts for SHE 
and OHE. Equations p3)) and (|14p enable one to separate 
the Hall effects, which lift the degeneracy of states with 
opposite helicities a — ±1 and vorticities I = ± and 
the diffraction effects described by Eq. (I15p . While the 
SHE turns out to be diffraction- independent, the OHE 
is crucially affected by the beam deformations upon the 
diffraction in a gradient medium. Due to this, calcula- 
tions of the OHE require numerical solution of the diffrac- 



tion Eq. (jlSp in each particular problem. 
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